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1 Introduction 

The origin of the large-scale structures of our universe can be explained by the presence of 
primordial density perturbations in the early universe, which were of super-horizon scales and 
almost scale-invariant. However, understanding how such density perturbations were seeded 
has remained a big mystery for cosmology. On the other hand, cosmic inflation [1, 2, 3] has 
been known to be able to generate super-horizon and (nearly) scale-invariant fluctuations for 
scalar flelds. (Early universe models that can generate such fleld fluctuations without having 
an inflationary era have also been studied, see e.g. [4, 5, 6].) Then now, one may rephrase 



the question about the origin of the large-scale structures as how such field fiuctuations were 
converted into the primordial density perturbations of the universe, which eventually lead to 
structure formation. 

An attractive mechanism for creating the primordial density perturbations from the field 
fiuctuations is the curvaton mechanism [7, 8, 9, 10], in which the curvaton field possessing 
large-scale field fiuctuations generates the density perturbations while it oscillates about its 
potential minimum and comes close to dominating the universe. The idea of the curvaton 
mechanism is thus simple, and also has the merit that when embedded in infiationary cosmol- 
ogy, it frees the infiaton from being responsible for generating the perturbations and therefore 
drastically relaxes constraints on infiationary model building. However, while there have been 
extensive works on the curvaton mechanism, most of the literature has focused on curvatons 
possessing rather trivial potentials (e.g. quadratic potentials), thus lead to common predic- 
tions for the density perturbations, such as large/small non-Gaussianity in the perturbations 
depending on whether the curvaton is subdominant/dominant upon decay. 

In this paper, we explore density perturbations sourced by a curvaton rolling along an 
arbitrary energy potential. The key feature of a curvaton potential which deviates from a 
quadratic one is that the curvaton field with large-scale field fiuctuations starts its oscillation 
at different times at different patches of the universe. Such behaviour contributes to the 
curvaton energy density perturbations in addition to that sourced directly from the original 
field fiuctuations. We find that such non-trivial conversion processes of the field fiuctuations 
into the density perturbations can lead to strong enhancement /suppression of the density 
perturbations of the universe, as well as for their non-Gaussianities. In particular, we will 
show that the non-Gaussianity parameter /nl can become large with either sign no matter if 
the curvaton dominates or subdominates the universe when it decays. Several previous works 
have considered specific cases, e.g., curvaton potentials with a non-quadratic polynomial 
term were studied in [11, 12, 13], and cosine type potentials in [14] which clearly pointed 
out that non-uniform onset of curvaton oscillations can have important effects on the density 
perturbations. Related discussions can be found in [15], though such effects were not explicitly 
taken into account. In this paper we carry out a systematic study of density perturbations 
from curvatons with generic potentials through analytic analyses. 

It is worthwhile to investigate curvaton potentials which deviate from simple quadratic 
ones, both from observational and theoretical reasons: Considering a curvaton whose field fiuc- 
tuations were generated during the infiationary era (which is the case studied in this paper), 
in order for the curvaton to produce the density perturbation spectrum which is red-tilted 
(i.e. with negative spectral index ng—l) as suggested by latest CMB observations [16, 17, 18] 
without relying on specific infiation mechanisms,^ the curvaton needs to be located along a 
potential with negative curvature during infiation (cf. (2.45)). This obviously suggests the cur- 



*For a non-tachyonic curvaton (i.e. a curvaton located along a potential with non-negative curvature during 
inflation), the spectral index of the resulting density perturbation spectrum is set by the time variation of 
the Hubble parameter during inflation as n^ — 1 > 2H/H^ (cf. (2.45)). Hence, assuming single-field canonical 



vaton potential to take non-trivial forms. Furthermore, explicit curvaton models constructed 
in the framework of microscopic physics can naturally possess intricate energy potentials, e.g., 
curvatons in supersymmetric models [20], a stringy curvaton model in [21] from a D-brane 
moving along the internal compactified space giving rise to a multi-dimensional periodic cur- 
vaton potential. [22] also discusses a string realization of the curvaton, where the potential is 
a linear combination of exponential terms. 

We also apply our general results to the case where the curvaton is a pseudo-Nambu- 
Goldstone (NG) boson [23] of a broken U(l) symmetry. The periodic curvaton potential 
necessarily possesses not only minima but also maxima, around which the potential curvature 
is negative so that the resulting density perturbation spectrum can become red-tilted. We 
explore the parameter space of NG curvatons producing a red-tilted density perturbation 
spectrum as suggested by observations, and show that the allowed regions for the inflationary 
energy scale and the reheating temperature (i.e. when the inflaton decays) are strongly 
constrained, unless the NG curvaton is located close to the maximum of its potential during 
inflation. In such hilltop cases, the resulting density perturbations are extremely enhanced, 
thus relaxing the bounds on the inflation/reheating scales. Non-Gaussianity of the density 
perturbations also increases in the hilltop limit to values 10 < /nl ^ 30 in most part of the 
allowed parameter window, which will be tested by upcoming CMB measurements. Such mild 
increase of non-Gaussianity towards the hilltop limit is actually a rather generic feature of 
hilltop curvatons, as will be shown in the paper. 

This paper is organized as follows. We derive analytic expressions for density perturbations 
from curvatons with generic potentials in Section 2. Then in Section 3, we move on to discuss 
special cases where the curvaton starts its oscillation from a flat region of its potential, 
including hilltop curvatons. It will be shown that in such cases, effects due to fluctuations 
of the onset of the oscillation become crucial. As a simple example of a curvaton model 
generating density perturbations that match with observational constraints, we investigate 
the case where the curvaton is a pseudo-NG boson in Section 4. We present our conclusions 
in Section 5. Discussions on dynamics of scalar fields in an expanding universe are provided in 
Appendix A. While we focus on sinusoidal curvaton oscillations in the main body of the paper, 
in Appendix B we further generalize our discussions to incorporate non-sinusoidal curvaton 
oscillations. Curvatons along linear potentials are also discussed in Appendix B. 



slow-roll inflation, the Lytli bound [19] relates the spectral index with the inflaton field (p range as 
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where Mp is the reduced Planck mass and Af the e-folding number. This shows that unless the curvaton is 
tachyonic during inflation, the WMAP [16] central value Ug ~ 0.96 requires a super-Planckian field range for 
the infiaton. 



2 Density Perturbations from Curvatons 

In this paper we focus on curvatons [7, 8, 9, 10, 24] whose field fiuctuations are generated in 
the infiationary era (though many features and results can be shared with curvatons with non- 
infiationary fiuctuation generating scenarios). Then a curvaton should be a light field so that 
it acquires field fiuctuations that are nearly scale-invariant and Gaussian. Its energy density 
is considered to be initially negligible, and after infiation the curvaton starts oscillating about 
its potential minimum and behaves as nonrelativistic matter. As its energy density grows 
relative to radiation sourced from the infiaton decay, the curvaton fiuctuations increasingly 
contribute to the density perturbations until the curvaton decays or dominates the universe. 
One of the main points of this paper is that the original field fiuctuations of the curva- 
ton 0" during infiation can give rise to perturbations of the curvaton energy density through 
perturbing the time when the curvaton starts its oscillation, as well as by directly perturbing 
its energy density p^. Let us schematically write down the resulting density perturbations of 
the universe as 

C~ci C2— — , (2.1) 

Pa -n osc 

where i^osc is the Hubble parameter at the onset of the curvaton oscillation. We have added 
a minus sign to the second term in the right hand side to indicate that larger i^osc gives an 
earlier onset of the curvaton oscillation, which more quickly redshifts away the curvaton energy 
density. (2.1) illustrates that direct energy perturbations 6pa can be cancelled by effects from 
SHosc, suppressing the linear order density perturbations with possible enhancement of the 
non-Gaussianity. Moreover, when effects from 6 Hose become dominant over that from 6p„, the 
resulting density perturbations become very different from cases with the familiar quadratic 
curvaton potential. However we should also stress that (2.1) is rather schematic, and the 
actual density perturbation spectrum is expressed as an involved combination of the two 
effects, cf. (2.39). 

In the following subsections, we explicitly compute the density perturbations from curva- 
tons with generic energy potentials. First deriving an expression of the form similar to (2.1) 
using the 5A/'-formalism [25, 26, 27, 28], we then analyze the curvaton dynamics and obtain 
analytic expressions for the density perturbations. 

We suppose that the Hubble parameter H = a/a (with a being the scale factor of the 
universe, and an overdot denoting a time derivative) is nearly constant during inflation, and 
that at the end of inflation, the inflaton suddenly turns into matter, i.e. its energy density 
starts decreasing as p^ oc a~^, and then at reheating, suddenly decays into radiation, i.e. 
starts decreasing as pr oc a~^. Curvature perturbations sourced from the inflaton are ignored. 
Furthermore, upon carrying out analytic estimations, we make the following simplifying as- 
sumptions for the curvaton: 

1. The curvaton a starts a sinusoidal oscillation at some time after inflation, from then 
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Figure 1: Schematic of the curvaton potential under consideration: The curvaton rolls along 
its energy potential following an attractor solution (2.24), then it starts oscillating about the 
(local) minimum. The potential can take arbitrary forms, but once the curvaton starts the 
oscillation (from ctosc) its amplitude soon decays and the potential is well approximated by a 
quadratic one. 



we treat the curvaton as matter. (Hence we are considering cases where the curvaton 
potential is well approximated by quadratic close to its (local) minimum, and that the 
curvaton oscillation happens mostly in the quadratic regime, cf. Figure 1. However, 
one can straightforwardly generalize the calculations to other cases, as is discussed in 
Appendix B.) 

2. The energy density of the curvaton until it starts oscillating is so tiny compared to the 
total energy density that its effect on the expansion of the universe is negligible at least 
until the onset of the curvaton oscillation. 

The curvaton acting as matter is also assumed to suddenly decay into radiation when H = Tfj, 
where we take the curvaton decay rate Fg. to be a uniform constant independent of the curvaton 
field value. 



2.1 Computing 5N 

The curvature perturbations C, from a curvaton can be computed using the 5A/'-formalism [25, 
26, 27, 28], as the difference in the e-folding number A/" at a constant energy hypersurface, 
among different patches of the universe which took different field values of the curvaton at 
the time when the CMB scale exited the horizon during inflation, i.e.. 
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Here, we are assuming a homogeneous and isotropic background, and values at the horizon 
exit are denoted by the subscript *. Expanding ({x) as 

Ck = j d'xe-'^-^a^), (2.3) 

their two- and three-point correlation functions take the form 

(a.a.) = {27:f5^'\kr + fc2)n(^i)> (2-4) 

(a^a^Cfca) = {2Tif6'^'\k^ + k^ + k^)B^{hM.h), (2.5) 

where ki = |fcj|. We define the power spectrum of the density perturbations P^ as 

Pdk) = ^n(^), (2.6) 

and also the non-linearity parameter /nl [29] as 

B^{h, k2, ks) = Ihdh, k2, ks) [P^{k,)P^{k2) + Pc(A:2)Pc(^3) + Pdks)Pc{k^)] , (2.7) 

which parametrizes the non-Gaussianity of the density perturbations encoded in the three- 
point function. From (2.2), these values can be evaluated in terms of the derivatives of A/", at 
the leading order in 6a^:k, as 

n-(^)V... (2.8) 

where P^o-, denotes the power spectrum of the curvaton field fluctuations defined similarly as 
in (2.6). Moreover, we have considered the intrinsic non-Gaussianity of the field fluctuations 
to be negligibly small and thus ignored the three-point function for Sa^,. 

2.1.1 Case trch < W 

Let us begin by studying the case where reheating (= inflaton decay, at treh) happens before the 
onset of the curvaton oscillation at tosc- In other words, the universe is radiation dominated 
when the curvaton starts oscillating. (Hereafter, values at time t, {i = reh, osc, etc.) are 
denoted with the same subscript.) We would like to derive the cr*-dependence of the e-folding 
number 

A^ = A4+A4, (2.10) 

where 

/■tosc rtdcc 

Ma= H{t')dt', K= H{t')dt'. (2.11) 



Here tdec denotes the time when the curvaton decays, i.e. when H becomes equal to V^ (which 
is independent of a^). After t^ec; the universe is filled with a single component, i.e. radiation, 
hence no extra 5M is produced. 

The energy density of radiation redshifts away as pr = —4:Hpr, thus 

A/; = /" "' H{t')dt' + [ °" H{t')dt' = const. + ^ In ^^^, (2.12) 
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where we have especially used Assumption 2. Since SM^H^ 
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Here we have introduced 
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Therefore, from (2.13) and (2.16) we arrive at 
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(where Mp is the reduced 
(2.13) 

(2.14) 

(2.15) 
(2.16) 

(2.17) 
(2.18) 



which reproduces the schematic expression of (2.1). We note that the second term in the 
right hand side that is due to the non-uniform onset of the curvaton oscillations, can become 
important for non-quadratic curvaton potentials.^ 

Computation of the second derivative of A/" is also straightforward. Since 

(2.19) 
(2.20) 
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''Strictly speaking, the first term also contains effects due to non- uniform onset of the oscillations through 
dcjosc/daif, as we will see later. Effects from 6Hosc displayed in (2.18) have not been taken into account 
explicitly in previous works giving general discussions on density perturbations in the curvaton scenario. 



2.1.2 Case treh > W 

The case of treh > ^osc can also be computed in a similar fashion as above. Since the curvaton 
energy density is assumed to be negligibly small compared to the total energy at the onset of 
the oscillation, the curvaton energy continues to be negligible until reheating. 
The linear order perturbation is 

where the coefficient of Inif^g^ differs from (2.18) due to different expansions of the universe 
(i.e. radiation or matter dominated) when the curvaton starts to oscillate. 
(2.20) remains the same, hence the second derivative of A/" is 

■a^ = WTWrW.) ■'its;*;! ('""■"- -'"^"=)- <"■"=" 

2.2 Curvaton Field Value at the Onset of Oscillation 

In order to evaluate the expressions derived above, we need to know the oscillation amplitude 
of the curvaton at a certain time, e.g., the curvaton field value at the onset of oscillation cXosc- 
Considering a curvaton rolling along an energy potential V which is a function only of a 
(hence we do not consider time-dependent potentials, though extending the discussion to 
such potentials is straightforward), then the curvaton dynamics before it starts its oscillation 
can be tracked by the attractor solution 

cH& ~ -1/', (2.24) 

where a prime denotes a derivative with respect to a, and c is a constant taking the following 
values depending on the expansion of the universe: 

3 (de Sitter) 
9/2 (matter domination) (2.25) 

5 (radiation domination) 

Detailed discussions on this approximation are provided in Appendix A. Here we assume 
that the approximation (2.24) holds until tosc, when the curvaton suddenly starts a sinusoidal 
oscillation, cf. Figure 1. Validity of this assumption will be discussed in the next subsection. 

Integrating 

da dt , . 

from the time of horizon exit until the onset of the oscillation, one obtains 

'"°=^ da ["°'^ dH 

— = const. - / ^, 2.27 



where c in the right hand side is the one right before tosc, taking either 5 or 9/2 depending 
on whether tj-eh < ^osc or treh > ^osc, respectively.^ We have made use of Assumption 2 upon 
ignoring the cr*-dependence of the other terms denoted as const. 

Supposing that ifosc is determined merely by o"osc, i-e. Hose = -f^osc(crosc) (cf. (2.33)), then 
after differentiating both sides of (2.27) by a*, one finds 



dan.r (^ 1 V'{aosc)dHosc\ ^ V'{aosc) ,^ ^^n 



da^ I c(c-3) H^^^ daosc j V'{a.. 

The value of ctosc itself can, in principle, also be computed from (2.27). Taking into 
account that H ~ ifinf = const, and c = 3 during inflation, and that H/H"^ = 3 — c for 
matter/radiation domination, one obtains for t^eh < tosc'- 

^ t 4.|.a-^V4f^-^V (2.2.) 

and for trch > ^o 



V 3HI, 27\Hl, HIJ 20\Hl^ Hi, 



'" ^- +4177^-7^1. (2.30) 



V 3/7^2 ' 97 \ t72 TT2 

where A/"* is the number of e-foldings during inflation between the horizon exit of the CMB 
scale and the end of inflation. Hence for i/^f ^ H^^^ (and further H^^^ ^ H^^^ for treh < W), 
we get 

^" ^* ' (2.31) 



V 3Hl, 2c{c-miJ 

where c is the value right before tosc- The two terms in the right hand side can be comparable, 
e.g., when setting the curvaton (effective) mass squared as rn? ~ ifj^f/lOO so that the spectral 
index becomes jn^ — 1| ~ 10"^, then Hl^^ r^ m"^ r^ H^^f/100 ~ Hf^JAf^:. Since H^^^ can be 
given by the curvaton potential as we will soon see, one can derive (Tosc as a function of cr* by 
solving (2.31). 

2.3 When the Curvaton Starts Oscillating 

The last piece of information needed to evaluate the density perturbation spectrum is when the 
curvaton starts its oscillation, i.e., -ffosc- From the necessary condition (A. 5) for the attractor 
solution (2.24) with (A. 6) to hold, one may guess that the oscillation starts when H'^ = \V"\/c. 
However actually, the curvaton velocity does not increase suddenly at H"^ = \V"\/c, but the 
curvaton stays on the approximation of the form H& oc —V for some more time. (See also 
later discussions in Subsection 2.5.) Moreover, the curvaton being on the attractor (2.24) 
does not necessarily imply that the curvaton oscillation has not yet started. 



^Strictly speaking, here we are further assuming that c stays constant around tosc, so that c is independent 
of the curvaton field fluctuations. So, e.g., the time of reheating should not be too close to tosc- 



Here, instead, we define the onset of tlie oscillation as when the time scale of the curvaton 
rolling becomes comparable to the Hubble time, i.e. 
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1, (2.32) 



where we have assumed that the (local) minimum of the curvaton potential is set to a = 0. 
Here it should be noted that modifying (2.32) by constant factors (e.g., choosing the right 
hand side to be, say, 2) in many cases results merely in modifying i^osc by constant factors, 
whose contributions to the density perturbations through 51nifosc drop out, cf. (2.18). Hence 
throughout this paper we fix the right hand side of the condition (2.32) to be unity. 

Then, supposing that the approximation (2.24) holds until tosc (which is a valid assumption 
in many cases since, as mentioned above, approximation of the form (2.24) does not break 
down suddenly at H^ = \V"\/c), one can rewrite (2.32) as 



H, 
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^'^"^"^^V (2.33) 



can 



Here c is the value right before tosc aiid we have removed the absolute value sign in the far right 
hand side by assuming that the curvaton potential is monotonically increasing (decreasing) for 
a > (<)0 (otherwise the curvaton following the attractor (2.24) stops rolling before reaching 
the origin). Thus we obtain 



1 aifosc 1 fV"(an.r) 1 



Hose da^sc 2 V V'(an<,r) cr„ 



(2.34) 



Together with (2.28) this yields 



dansr ( 1 f croscVjcTosc) .\\ ^ V^Ooscl ^2 35^ 



da, { 2(c-3) V V'iaosc) J } V'{a,) 

2.4 Results 

Putting together the above results, and also adopting 

P<7osc = V{aosc), (2.36) 

we obtain the final expressions for the density perturbations. Let us define 

where c = 5 for trch < W and c = 9/2 for t^^h > ^osc- Note that X denotes effects due to 
fluctuations of the time when the curvaton starts its oscillation 6Hn^r. 
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Then using Vg^ = if*/27r, the hnear order density perturbation spectrum is 
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The second order perturbations in terms of the non-hnearity parameter /nl are 
5 9W fdM^~^ 



(2.38) 
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Thus we have managed to write down curvature perturbations from a curvaton in terms of 
(^"5 CTosc, cr*), where Oo^c is determined by a^, through (2.31). 

The curvaton energy fraction upon decay r (2.17) can be estimated by assuming that the 
universe suddenly switches between matter /radiation domination at reheating (which happens 
when H = F,^, with T^ the inflaton decay rate) and also when the curvaton dominates the 
universe (if it ever does) , 
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Note that we have again especially made use of Assumption 2. The first and second terms in 
the Max. parentheses correspond to the curvaton being subdominant and dominant at decay, 
respectively. Also, the Min. parentheses are due to whether the onset of oscillation is after or 
before reheating. 
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example: Quadratic Potential 

As a simple example, let us consider the familiar case of a quadratic potential V cc a"^. One 
sees that X vanishes, i.e. the onset of the curvaton oscillation is uniform, and 

dAf 2r 1 5 / 4 8 \ , , , 

/nl = T7T -3 + - + ^^r ' (2-42) 



da^ 4 + 3r(j,' ' 12 V r 4 + 3r 

where the well studied results are reproduced. Here, in particular, large non-Gaussianity 
/nl ^ 1 is sourced only by curvatons decaying while r ^ 1. 

2.5 Validity of the Analytic Expressions 

Upon deriving the above analytic results, we have adopted simplifying assumptions such as 
sudden decays, and also that the curvaton suddenly starts a sinusoidal oscillation when (2.32). 
Here we comment on the validity of such assumptions and the analytic results. 

First we should remark that upon obtaining the expression (2.33) for Hose, we have as- 
sumed the curvaton to follow the attractor solution (2.24) until the onset of the oscilla- 
tion (2.32). As is discussed in the appendix, this procedure is vahdated under the necessary 
condition (A. 5). Thus one should be careful when dealing with curvaton potentials along 
which (A. 5) breaks down significantly before the oscillation. There may also be cases where 
(A. 5) once breaks down but recovers afterwards, such that the curvaton once deviates from 
the approximation (2.24) but comes back again later on, if, say, the curvaton potential pos- 
sesses a plateau region. However, we also note that the breakdown of (A. 5) before the onset of 
the oscillation does not necessary imply the breakdown of the analytic results derived above. 
When \V" /H'^l becomes of order unity, the parameter c in the approximation (2.24) no longer 
takes constant values (2.25), but becomes time-dependent as can be seen from (A. 4). Hence 
using the time-dependent solution of (A. 4) for c (given that the solution exists^) may give 
more accurate estimations, but as long as the time scale of the variation of c is of order the 
Hubble time, our discussions with constant c should not be drastically modified. Moreover, 
unless (A. 5) breaks down long before the onset of the oscillations, the deviation of the cur- 
vaton dynamics from (2.24) with constant c is mild and the above analytic results stay valid. 
We will see in Subsection 4.1.2 an explicit example where the condition (A. 5) clearly breaks 
down by the time when the oscillation starts, but still the above analytic results giving good 
estimations of the density perturbations. 

Upon estimating Po-osc, we have adopted (2.36) and ignored the curvaton's kinetic energy. 
This procedure is validated in many cases where the potential energy is larger than the kinetic 
energy at the onset of the oscillations. Moreover, as long as the ratio between the potential 
and kinetic energies is a constant, considerations of the kinetic energy drop out through 



^Note that (A. 4) admits solutions of c for largely negative V" /H^, but has no solutions for largely posi- 
tive V"/H^. 
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(51npo-osc- However, there may be cases for which it becomes important to take into account 
the curvaton's kinetic energy at tosc- 

It was also assumed that even if the curvaton potential deviates from a quadratic one 
(giving rise to 6 Hose), once the oscillation starts its amplitude soon decays and the curvaton 
starts a sinusoidal oscillation about its minimum where the potential is well approximated 
by a quadratic one. However, we should also note that if the curvaton oscillation starts 
with an amplitude large enough so that its potential deviates from quadratic for a while 
after tosc, then the curvaton energy density does not behave as matter but decays as p^ oc a~" 
with n 7^ 3. This can modify our analytic expressions, in addition to changing r and thus 
the amplitude of the perturbation spectrum. In cases where the oscillation history is rather 
trivial (e.g. V = ^m?(j'^ + Ao"^ giving n = 4 — ;■ n = 3), extending the discussions above is 
rather straightforward, as is shown in Appendix B. However for more complicated cases (e.g. 
curvaton oscillation in a potential with superimposed periodic modulations, which can give 
rise to a time- dependent n), some parts of our results may require major corrections. We 
leave this for future work. 

2.6 Spectral Index 

Although the perturbation amplitude crucially depends on the curvaton dynamics after in- 
flation, the spectral index 

is determined only by the information at the horizon exit during inflation. Note that the 
linear order density perturbation is expressed by 

TT— = (function of r, cTosc) x — - — -, (2.44) 

where r and o"osc are independent of the comoving wave number k. Hence 

d f H, Y 2V"ia,) ti. 



where we have made use of d\n.k ~ H^dt and 3H^:&^: ~ — V(a^,). One clearly sees that for 
inflation with H^, = const., the observationally suggested red-tilted perturbation spectrum 
requires the curvaton to be tachyonic during inflation. See also discussions in Footnote 4. 

3 Curvatons Along Flat Potentials 

The results (2.39) and (2.40) in the previous section show that additional contributions to the 
density perturbations from ^i^osc can suppress/enhance the perturbation amplitude as well 
as the non-Gaussianity. As an interesting case where the SHosc contributions to the density 
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perturbations become important, in this section we focus on a curvaton starting its oscillation 
from a flat region where the potential tilt is tiny. 

Given that the first derivative of the potential at the onset of oscillation is small to enough 
to satisfy 



a„s.r'(a, 



osc ^ \^ osci 



T^'(^osc) 

then one finds 



^"(aosc)V^(aosc) 



1 



^'(^osc)^ 



O'osr-V^"(crn 



V"(a_)2 



V'{a^,,)V"'{a, 



(^oscVicT*) 



Via,, 



>1, 



2(c - 3) V\a, 



> 1, X\a, 



1 CFa^cVian 



2(c-3) y'(aosc)2 



(3.1) 



(3.2) 



and the equations (2.39) and (2.40) are approximated by 






3r V'[ 



cr„ 



/i 



NL 



4 + 3r(ToscV^'(a.)' 
5(4 + 3r)aoscr'(cr. 



(3.3) 



(3.4) 



18r V'(aosc) 

Note that here, |/nl| ^ 1 is guaranteed (even when r ^ 1). Furthermore, when V{a) 
is a monotonically increasing function around its origin, i.e. V > (<)0 for a > (<)0, then 
/nl oc (negative factor) x V"{a^). In other words, with constant H inflation, a red(blue)-tilted 
perturbation spectrum is always accompanied by a largely positive (negative) /nl- 

We should remark that under the condition \aV"/V'\ ^ 1 in (3.1), (A. 5) no longer holds 
at tosc- This may lead to some deviation of the curvaton dynamics from the approxima- 
tion (2.24), giving rise to errors in the estimations. However, as we will see in Section 4 when 
we study explicit cases, the overall behaviour of the resulting density perturbations discussed 
here are still valid as long as the breakdown of the approximation (2.24) is mild. (See also 
discussions in Subsection 2.5). 



3.1 Hilltop Curvatons 

Let us further focus on a curvaton located at the hilltop, i.e. a curvaton whose potential is 
well approximated by 



V{a) = Vo — -rn^icr — (Jq 



(3.5) 



around (Tosc and a*, where m, ctq, and Vo(> 0) are constants. Such case will be important when 
we study pseudo-Nambu-Goldstone curvatons in the next section. Without loss of generality, 
hereafter we assume < o-qsc < ct* < o"o- 

Now the flatness conditions in (3.1) translate into 



cr„ 



>1, 



Vo 



>1, 



(^0 



0"n 



m" (Tn 



CTOJ 



(3.6) 
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under which (3.3) and (3.4) are 



V, 



1/2 



3r ctq 



o-n 



H: 



inf 



/l 



NL 



4 + 3r (To - (J* 27rcro 

^ 5(4 + 3r) gpsc 

18r ctq - o-Qsc ' 



with spectral index 



2 m^ 



(3.7) 
(3.8) 

(3.9) 



Here we note that throughout this subsection, inflation with constant H is assumed. 
The Hubble parameter at the onset of oscillation (2.33) is 



m 



m^{oQ - o-Q 



C(Tn 



< Hini^ 



(3.10) 



Thus (T* and (Jn^r can be 



where the inequality follows from (3.6) and In^ — 1| = (9(10 

related by solving (2.31), which in the hilltop limit, i.e. a* — )■ a^ (from below), gives 

(7q — (7* \ 1 



In 



o"o - o-Q 



CTn 



2(c-3)ao 



O-n 



(3.11) 



where c = 5 for trch < i^osc and c = 9/2 for treh > ^osc- Since the left hand side is logarithmic, 
one sees that o"osc approaches ctq much slower than a^. does, cf. Figure 2. Thus, in the hilltop 
limit, P( (3.7) blows up due to the enhancement factor (ctq — crosc)/(o"o — a*), while /nl 
(3.8) increases slowly. One can see from this example that the evaluation of o"osc is crucial for 
correctly estimating the density perturbations. The extreme amplification of the perturbation 
amplitude can be understood as curvatons taking longer time to start its oscillation as one 
approaches the hilltop limit. 




Figure 2: cXosc in terms of a^, obtained from (3.11) with c = 5 (red) and c = 9/2 (blue) 
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One may worry, in the hilltop limit, that the curvaton during inflation jumps over the 
hill due to quantum fluctuations Act = ±iJinf/27r, and domain walls be generated. Here we 
note that if the hilltop curvaton generates density perturbations that satisfy the observational 
constraints from the WMAP 7-year data [16] Vc^ ~ 2.4 x 10~^ and /nl < 74 (for local-type 
non-Gaussianity which is important for us here), the condition for the absence of domain 
walls is satisfled, 

^(ao - a.)-' ^ ^Vl^'fNL < 0.004. (3.12) 

In 5 ^ 

We also have to check whether quantum fluctuations lead to a random walk of the curva- 
ton. Further using the WMAP7 central value Ug ~ 0.96, one sees that the classical rolling 
dominates over the quantum fluctuations, 

< 0.2. (3.13) 



H^ 



inf 

~27 



cr^ 



m 



inf 



Hi, ^ 12 P^^^7n l 
2ti V'ia^ ~ 5 l-n. 



Since the onset of the curvaton oscillation is delayed as one approaches the hilltop limit, 
it becomes easier for the curvaton to dominate the universe before it decays. However, we 
should also note that the curvaton energy fraction when it starts to oscillate, 

"- - "■ '' ^"^ """ St^J;Wx™00. (3.14) 



3M|ffi, 4 + 3r 5 1 - n, 3Af ^i/? , ~ 4 + 3r SM^H^ 

can become as large as unity. In such case, the above estimations can break down since upon 
deriving the analytical expressions, we have ignored the effects of the curvaton energy density 
on the evolution of the universe before the curvaton starts its oscillation, cf. Assumption 2. 

4 Case Study: Pseudo-Nambu-Goldstone Curvatons 

As a simple example of a curvaton model, let us examine the case where the curvaton is a 
pseudo-Nambu-Goldstone (NO) boson of a broken U(l) symmetry. The approximate sym- 
metry suppresses the curvaton mass, while the periodicity of the U(l) provides both (local) 
minima and maxima along the curvaton potential, thus allowing the curvaton to produce a 
blue-tilted as well as a red-tilted density perturbation spectrum. 
In particular, we study the NG curvaton potential of the form 



V{a) = A' 



a 

cos: 



(4.1) 



where / and A are mass scales. Without loss of generality, we set the fleld value of the 

curvaton at horizon exit during inflation to lie within the range < o"* < vr/, and consider 

its oscillation about the origin a = 0. Supposing that the coupling of the NG curvaton with 

its decay product is suppressed by the scale of symmetry breaking /, then the decay rate of 

the curvaton takes the value 

(3 m^ 13 K^ 



167r P 167r p 
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(4.2) 



where m is the mass at the miniinuiii, i.e. m^ = V"{0), and /3 is a constant whose precise 
value depends on the decay modes. Throughout this section we fix /3 to unity. (Cases with 
smaller /3 will be discussed towards the end of Subsection 4.2.2). Furthermore, we consider 
infiation with constant H, and assume the pivot scale to have exited the horizon 50 e-foldings 
before the end of infiation. 

We will see that an NG curvaton away from the hilltop requires infiation and the reheating 
scales to lie within a narrow range in order to generate density perturbations with spectral 
index of order n^ — 1 ~ —10^^. However, we also show that as one approaches the hilltop 
limit (T=K — )■ vr/, the allowed infiation/reheating scales extremely broaden and that the non- 
Gaussianity lies in the range 10 < /nl ^ 30. 

For an investigation of NG curvatons located close to their potential minimum, see also [23] . 

4.1 Density Perturbations from NG Curvatons 

Let us first investigate density perturbations generated by NG curvatons, and compare the 
analytic estimations in the previous sections with numerically computed results. We start 
from NG curvatons located away from the hilltop of the potential, then move on to the 
hilltop limit. 

4.1.1 Non-Hilltop Region 

Firstly, we choose the potential parameters and the infiation/reheating scales such that the 
curvaton with cr* = ^Tcf generates a density perturbation spectrum whose amplitude takes 
the COBE normalization value 

V^ ^ 2.4 X 10"^ (4.3) 

with the spectral index lying at the central value of the WMAP7 [16] bound (for a power-law 
spectrum with no tensor modes) 

Us ^ 0.96. (4.4) 

Then, considering the case where tosc < ^reh, and that the curvaton dominates the universe 
before it decays, we adopt the following parameter set (which lies at around the center of the 
allowed window in Figure 11 in the next subsection): Hmf = 10^^'^ GeV, p/^^ = 10^^ GeV (i.e. 
r^ ^ 9.73 X 10-^Mp), / ^ 3.36 x lO'^Mp, and A ^ 3.56 x lO'^Mp. 

In order to analytically estimate the density perturbations, we need to solve (2.31) for 
obtaining cXosc- Introducing 



a 



one obtains 



In 



tan (aosc7i"/2) 



tan {a^7r/2) 



a = — (4.5) 



A4 A"^ 1 ttoscvr ,. „^ 
(4.6) 



?,Hl,P 2(c-3)sin(ao,,7r)' 
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whose solution for A/"* = 50 and c = 9/2 in the region 0.01 < a* < 0.99 is approximated by 
(cf. Figure 3) 

aosc ~ 0.195a* - 0.0993a^ + 0.846a^ - 1.34a^ + 1.26a^ (4.7) 

Then by plugging (2.41), (4.1), and (4.7) into (2.39) and (2.40), one obtains the linear order 
density perturbations and its non-Gaussianity as functions of a* = o^jixf. 

0.8 
0.6 
0.4- 

0.2^ 

Figure 3: o"osc in terms of a* for the NG curvaton, cf. (4.7). 



We have also numerically computed the density perturbations by solving the equation of 
motion of the curvaton. After 50 e-foldings of inflation with constant H since the exit of the 
pivot scale, the infiaton energy is transferred to infiaton matter decaying as pm oc a~^, which 
eventually decays into radiation (p^ oc a~^). Sudden decays are adopted for the infiaton and 
the curvaton, i.e. they instantly decay into radiation when H = T. 

We run the numerical computations with different a* (the curvaton velocity at the horizon 
exit is set to the slow-roll attractor value, which is uniquely determined by cr*), and evaluate 
the number of e-foldings obtained by the time the curvaton decays. Thus evaluating the 
differences in e-folding numbers AN' due to shifts in the initial condition Act, we numerically 
obtain dM /da^,. Second order perturbations as well as the spectral index are computed in 
similar manners. 

In Figures 4 - 7 we lay out the results in terms of a*, where the solid lines are the analytic 
estimations and the dots are the numerical results. Blue lines and dots correspond to our case 
with Fg- = -^yi where the WMAP values (4.3) and (4.4) are realized at a^ = 3/4, and for 
comparison we also show cases with Fo- = 10^ ^ li^i^Ta" (red) and T„ = 10^"^ ^ jki^j^ (green). 
(Note that Figure 4 has only one line since the spectral index is independent of the decay 
rate.) One sees that the results asymptote to that of the familiar quadratic curvaton (2.42) 
as a* — ;■ 0. On the other hand, the linear order perturbation as well as the non-Gaussianity 



increases as one approaches a* — )■ 1, smoothly connecting to the hilltop limit which we will 
soon study. The curvaton energy fraction at decay r increases with a^,, and the density 
perturbations become r-independent for r ^ 1, as can be seen from the overlap of the blue 
and green lines towards a* — )■ 1. 




1.0 tt/ 




Figure 4: Spectral index. 



Figure 5: Energy fraction at decay. 




1000- 




Figure 6: Power spectrum. 



Figure 7: Non-Gaussianity. (The green line 
goes negative in the covered region.) 



4.1.2 Hilltop Region 

As we have seen in Subsection 3.1, an extreme enhancement of the density perturbations 
occurs in the hilltop region. Let us now choose parameters such that density perturbations 
with (4.3) and (4.4) are realized at a* = (1 — 10~^)7if. Considering the case where tosc < ^rch 
and the curvaton dominates the universe before it decays, we adopt the following parameter 
set (which lies at around the center of the allowed window in Figure 15): //mf = 10^ GeV, 
pj/^ = 109 GeV (i.e. F^ « 9.73 x IQ-^OMp), / « 1.57 x IQ-^Mp, and A « 1.23 x IQ-^Mp. 
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Introducing n^ as 

a, = 1-10-"*, (4.8) 

the solution of (4.6) in the region of our interest 2 < n* < 10 is approximated by 

ttosc ~ 1 - (0.416 + 0.0120n,)n;^-^^ (4.9) 

which shows that aosc approaches unity much slower than a^, does (cf. Figure 2). Then 
one can evaluate the density perturbation spectrum as was done for the non-hilltop region. 
The results are shown in Figures 8 - 10, where again the solid lines denote the analytic 
estimations and dots the numerically computed results. The curvaton well dominates the 
universe before it decays, i.e. r ::^ 1, in the displayed region 5 < n,, < 9. Here we should 
remark that, in the hilltop limit, the attractor (2.24) (mildly) breaks down before (2.32) is 
realized, hence (Josc (4.9) estimated from (2.33) contains error of order AcXosc/cosc = 0{1).^ 
This gives rise to errors of the similar order in the analytic estimations, which is clearly 
seen in Figure 10. Therefore, we have also calculated semi-analytic results by substituting 
numerically computed ctosc (defined by (2.32)) into the analytic expressions in Subsection 2.4. 
The semi-analytic results are shown in the figures as dashed lines, which match well with the 
numerical computations. (The dashed line is absent in Figure 8 since the spectral index is 
evaluated independently of o"osc)- 

The figures show that, as was discussed as generic features of hilltop curvatons in Sub- 
section 3.1, the linear order density perturbations blow up in the hilltop limit, while the 
non-linearity parameter /nl increases slowly. The spectral index, since it is determined by 
the potential curvature, asymptotes to a constant value as one approaches the hilltop. How- 
ever in the extreme hilltop limit the finite size field fluctuations 5a^ = ifinf/27r can no longer 
be treated as infinitesimal, giving rise to the tiny variations of the numerical results seen in 
Figure 8. 

4.2 Parameter Space 

Now that we have studied the behaviour of density perturbations from NG curvatons, let 
us explore the parameter space which allows generation of observationally consistent density 
perturbations. Out of the five free parameters (/, A, -ffmf, -ffreh, cr*), two are fixed from the 
COBE normalization (4.3) and the WMAP7 central value for the spectral index (4.4). We 
then seek constraints on the inflation/reheating scales under various a^,, and discuss cosmology 
with NG curvatons. 

^The value of —V' /Ha (which in the analytic estimations we have treated as a constant c — 9/2) becomes 
as large as ~ 30 right before the onset of the oscillation (2.32). (This can be estimated by solving the 
equation (A. 4) in terms of c.) Nevertheless, the error in CTosc is of order unity since —V'/Ha being large is 
limited to times close to iosc- 
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Figure 8: Spectral index. 
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Figure 9: Power spectrum. 



Figure 10: Non-Gaussianity. 



21 



4.2.1 Requirements 

Let us begin by laying out the requirements for a consistent curvaton scenario. Firstly, we 
assume the curvaton energy density to be neghgibly small during inflation, 

^("*) «1, (4.10) 



^M^M, 



and until the curvaton starts its oscillation (cf. Assumption 2) 



^M^Hl^ 



< 1. (4.11) 



Also, quantum fluctuations during inflation should not make the curvaton jump over its 
potential minimum in order to avoid the resulting density perturbations from being highly 
non-Gaussian, or over the maximum to avoid domain walls, 

^ « ^. « ^/ - ^- (4.12) 

We further require the curvaton's classical rolling to dominate over the quantum fluctuations 
during inflation (so that we can have definite predictions for the perturbation spectrum^), 

^ ^'"' <1, (4.13) 



271 Via, 



where the slow-roll approximation for the curvaton during inflation is verified by (4.10) and 
(4.19). The curvaton decay should happen after reheating 

r^ > r., (4.14) 

and also after the onset of the curvaton oscillation 

^osc > r^, (4.15) 

but at the latest before the Big Bang Nucleosynthesis (BBN) at temperature Tbbn ~ 1 MeV 
[30, 31, 32, 33], 

SM^Tl > {lMeV)\ (4.16) 

We also require the oscillating curvaton's mass to be larger than its decay temperature, in 
order to avoid possible backreaction effects to the curvaton's perturbative decay (see e.g. 
[34, 35, 36] for discussions on this issue). Assuming instant thermalization, this condition is 
written roughly as 

V"{0) > (SM^Tiy/^. (4.17) 



'The randomized case of NG curvatons are discussed in [23]. 

22 



The resulting density perturbations should satisfy the COBE normalization, 

V(^ ^ 2.4 X 10"^ (4.18) 

and for its spectral index we adopt the WMAP7 central value (4.4), which through (2.45) 
translates into (note that we are supposing inflation with H = const.) 

?^«-0.04. (4.19) 

We note that when (4.19) is satisfied at a* ~ 7r//2 with |cos(cr*//)| ^ 1, then a rather 
large A^//^ can give rise to a large running of the spectral index and may contradict with 
observations. However, when fixing n^ — 1 ~ — 10~^, the allowed parameter window becomes 
small anyway as cr* approaches the inflection point 7r//2 (as we will soon see), thus we merely 
require (4.19) for the spectral tilt. Bounds on local-type non-Gaussianity from the WMAP7 
(at 95% CL) are 

- 10 < /nl < 74. (4.20) 

Constraints on primordial gravitational waves set an upper bound on the inflationary 
energy scale. The 7-year WMAP+BAO+/7o gives Vt/V^ < 0.24 (95% CL), i.e., 

ffinf < 1.3x lO^^GeV. (4.21) 

However, we note that when the inflation scale is high enough to (come close to) saturate 
the bound (4.21), then depending on the inflationary mechanism, one can expect to have 
a contribution to the spectral index from a non- vanishing H/H^ (see also Footnote 4), as 
well as the central value of the spectral index bounds (4.4) being slightly shifted (to about 
0.97 with the 7-year WMAF+BAO+i^o)- Here, as we have stated at the beginning of this 
subsection, we adopt the value (4.4) and a H = const, inflation, and simply use (4.21) as an 
upper bound on the inflationary scale. 

Finally, reheating, i.e. inflaton decay, should happen subsequent to inflation, 

F^ < //inf. (4.22) 

4.2.2 Windows for Inflation/Reheating Scales 

Among the above constraints, (4.18) and (4.19) can be used to flx the parameters / and A. 
Here, note that ctosc// is estimated as a function of cr*// by solving (4.6) under AC = 50, 
setting c = 5 (or 9/2) for tosc > (or <) treh, and choosing A'^/ pH^^^ as a function of 0"=,,// from 
(4.19). 

Thus we obtain windows in the parameter space that satisfy all the above constraints 
(4.10) - (4.22) in four regimes, categorized by r ^ 1 and tosc ^ ^reh- The results are displayed 
in Figures 11-16, in the ifmf - p^eh plane under various values of a*//. The allowed windows 
are shown as colored regions, where the yellow (green) regions correspond to that for r > 1 
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with tosc > i^reh (^osc < ^reh), and the blue (red) regions for r < 1 with tosc > Ueh (tosc < treh), 
respectively. The right side edge of each box, i.e. the maximum displayed value of -ffinf, is 
set by the gravitational wave bound (4.21), also the gray shaded region is excluded since 
the inflaton should decay after inflation (4.22). Differently colored boundaries of the allowed 
regions denote different constraints being saturated, where purple dashed: curvaton being 
subdominant at the onset of oscillation (4.11), green: curvaton decay after reheating (4.14), 
orange dashed: BBN (4.16), and red: upper bound on /nl in (4.20). Furthermore, the blue 
lines denote where r = 1. The inequalities ^ and 3> in the above conditions are relaxed 
to < and > upon obtaining the allowed windows in the figures, hence in this sense the 
displayed constraints are conservative bounds. ^° Moreover, we note that we have examined 
the constraints independently in four distinct regions tosc ^ ^reh (^ tj-eh) and r ^ 1 (^ 1), 
and displayed (two of) them together in each figure. This is why the displayed boundaries bend 
sharply when crossing tosc = treh and moving to differently colored regions. However, we expect 
that such simplifications and also errors due to approximations used in the analytic analyses 
(including that upon estimating cTosc in the hilltop limit, see discussions below (4.9)) are minor 
for obtaining order-of-magnitude estimations of the inflation/reheating scales compatible with 
NG curvatons. We also show contour lines of the curvaton decay temperature for r > 1 cases, 
where the numbers in boxes denote (SM^F^)^/^ [GeV]. (Note that the contour lines are 
actually relevant only on the colored allowed regions.) 

First of all, requiring a red-tilted density perturbation spectrum, the NG curvaton at 
horizon exit should be located at 1/2 < a,, = o^^j-n] < 1. Constraints on NG curvatons away 
from the hilltop can be seen in Figures 11 and 12, where the representative value a* = 3/4 is 
taken. The dominant and subdominant cases share the same regions in the ifinf - Pj^f plane, 
bounded by (4.21) and (4.22). Here, one sees that NG curvatons away from the hilltop require 
rather high inflation/reheating scales. This is because the spectral index In^ — 1| ~ 10~^ 
indicates a curvaton's effective mass as large as |meffp ~ 10~^i7^f, which starts the curvaton 
oscillation soon after inflation ends.^^ Unless the inflation and reheating scales are high 
enough to provide a sufficiently long radiation-dominated era for the oscillating curvaton to 
dominate the universe, the curvaton cannot source substantial density perturbations. This 
is especially the case as one approaches the inflection point of the potential, i.e. a* -^ 1/2, 
since then | cos((j*//)| — )■ and the tilt of the potential increases, pushing i^osc (2.33) as well 
as the inflation/reheating scales even higher. 

On the other hand, as one approaches the hilltop limit a* — )■ 1, the onset of the oscillation 
is delayed and curvaton domination is allowed with lower inflation/reheating scales. Moreover, 



^°Even under the relaxed conditions, we especially note that ojJ^|^2 < r is guaranteed by the conditions 
(4.14), (4.15), and the analytic expression for r (2.41). 

^^Strictly speaking, as we have emphasized in the previous sections, it is not the potential's curvature but 
rather its tilt that determines the onset of the oscillation. However, away from the inflection point ct/tt/ = 1/2 
or the hilltop ct/tt/ — 1, the curvature and tilt are of the same scale, i.e. \V"\ ~ V /cr, hence either can be 
used for estimating i/osc- 
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the linear order density perturbations obtain an enhancement factor ( 1 — ttosc) / ( 1 ~ c^* ) ? which 
should be compensated by a smaller -f/^mf// (cf. (3.7)). This translates into a larger A through 
the fixed spectral index (4.19), thus further increases the curvaton energy density fraction, 
cf. (3.14). Figure 13 with a* = 1 — 10^^ shows that lower infiation/reheating scales are 
compatible with hilltop NG curvatons, thus broadening the allowed window. One also sees 
that now our requirement of the curvaton being subdominant until the onset of oscillation 
(4.11) becomes important, cutting off the right edge (i.e. high inflation scales). 

As we have noted in Subsections 3.1 and 4.1.2, non-Gaussianity also increases, though 
mildly, towards the hilltop. In the allowed region in Figure 13, the non-Gaussianity is 
/nl ~ 10, except for close to the left edge where r ~ 1 and hence larger /nl. (Here, upon 
evaluating the explicit value of /nl, we have carried out numerical calculations in order to 
avoid errors rising from the approximations used in the analytic estimations.) For a subdom- 
inant curvaton whose non-Gaussianity is further enhanced, the observational upper bound on 
/nl (4.20) becomes a severe constraint as is seen in Figure 14. Eventually the allowed window 
for subdominant curvatons vanishes as one approaches the hilltop and the non-Gaussianity 
increases. When log]^o(l ~ o;*) is about —8, the window exists only for dominant curvatons, 
which is shown in Figure 15. Here, the left edge is now the non-Gaussianity bound. 

Further approaching the tip, the BBN constraint (4.16) kicks in from low -ffinf, cf. Figure 16 
where «=„ = 1 — 10^^^. The non-Gaussianity in most of the allowed region of Figure 16 (except 
for close to the bottom edge) is /nl ~ 30. As one goes even closer to the hilltop, the BBN 
constraint goes towards higher H^^f while the constraint (4.11) goes towards lower ifjnf, thus 
removes the allowed window completely by logio(l — a*) ~ —12. 

When the curvaton decay rate is further suppressed, i.e. when /3 in (4.2) is way smaller 
than unity, the overall behaviour of the allowed window stays the same except that the 
curvaton decaying temperature is lowered correspondingly. For instance, when /3 = 10~^, 
the blue lines in the figures where r = 1 are shifted towards lower H^^f by about one order 
of magnitude (or towards lower p^^^ by about two orders). The red lines (/nl bound) shift 
similarly while the purple dashed lines (4.11) are independent of /3. Therefore the allowed 
windows for dominant curvatons slightly broaden in Figures 11, 13, 15, while the windows 
for subdominant curvatons simply shift towards lower i^inf, Pj^f . However, since the BBN 
constraint becomes severer for lower curvaton decay temperatures, the allowed window for 
dominant curvatons vanishes at logio(l — a*) ~ —10. 

Let us comment on implications of the present scenario for the origin of the baryon asym- 
metry and dark matter. Focusing on the case of the dominant curvaton, the adiabatic density 
perturbation becomes fixed when the curvaton dominates the universe. In order to avoid 
generating too large isocurvature perturbations, therefore, the baryon as well as dark matter 
abundance must be fixed after that. The thermal leptogenesis scenario [37] requires a rela- 
tively high curvaton decay temperature Tdcc ^ 10^ GeV. In the present scenario, T^ec exceeds 
10^ GeV for 1 — a^, > 10~^. As «=„ further approaches 1, the decay temperature decreases and 
other baryogenesis scenarios are needed. For instance, non-thermal leptogenesis [38] is known 
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to work for T^^c ^ 10^ GeV, which will be viable for 1 — a,, > 10^^. For the Affleck-Dine 
mechanism [39, 40] to work, the following condition must be met, -ffdom > "^adj where i?^dom 
denotes the Hubble parameter when the curvaton dominates the universe, and wiad(^ ^^3/2) 
denotes the mass of the Affleck-Dine field. The above condition can be met for various a* as 
long as the reheating scale is not so low. The electroweak baryogenesis requires the curvaton 
decay temperature to exceed of order 100 GeV, which is translated into 1 — a^ > 10~^. For 
1 — q;^, < 10~^, the decay of a may be able to generate the baryon asymmetry in a way similar 
to Ref. [41], using the R-parity violating operators. On the other hand, there are many dark 
matter candidates whose abundance is fixed after the curvaton dominates the universe. For 
instance, the (non-)thermal relic of the WIMPs and the QCD axion would fit with the present 
scenario. 

Summarizing, the allowed parameter window for NO curvatons in the i^inf - Pj^f plane 
broadens significantly in the hilltop regime, especially for 1 — lO"'^ < cr*/7r/ ^ 1 — 10~^^. In 
such hilltop regime, the NG curvaton mainly dominates the universe before it decays, but still 
produces non-Gaussianity in the density perturbations of size 10 < /nl ^ 30. We reiterate 
that this prediction on /nl is free from the errors due to approximations used for the analytic 
estimations, although the constraints displayed in the figures contain errors corresponding to 
order unity factors. 
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Figure 11: crjnf = 3/4, r > 1. 



Figure 12: ajnf = 3/4, r < 1. 
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Figure 13: a^/nf = 1 - 10-^ r > 1. 
/nl ~ 10 in most of the allowed region. 



Figure 14: a^/nf = 1 - 10-^ r < 1. 
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Figure 15: ajnf = 1 - 10-^ r > 1. Figure 16: ajnf = 1 - 10"", r > 1. 
/nl ~ 20 in most of the allowed region. /nl ~ 30 in most of the allowed region. 
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4.3 Discussion on Hilltop Curvatons 

Let us comment on initial conditions of hilltop curvaton models. So far we have simply 
assumed that the curvaton sits near the local maximum. This may be justified if the curvaton 
is charged under some symmetry, which is unbroken during early stage of the inflation or pre- 
infiationary epoch. This ameliorates the amount of fine-tuning of the initial condition for the 
hilltop curvaton. On the other hand, in the NG curvaton model, there is no special point in its 
field space since the curvaton has an approximate shift symmetry. Then, a priori, there is no 
particular reason to believe that the curvaton initially sits near the local maximum. However, 
we do not know whether the conventional naturalness argument can be really applied to 
the determination of parameters in the early universe. For example, the parameters may be 
determined in a similar fashion as in the anthropic argument for explaining the present value 
of the cosmological constant [42]. Suppose that the prior probability of the initial position 
of the NG curvaton has a fiat distribution in the field space. Assuming that all the other 
parameters are fixed, it is conceivable that only the universe with the curvaton initially sitting 
near the local maximum is anthropically selected, if it does not lead to the habitable universe 
otherwise. For instance, this is expected to be the case, if the infiaton mainly decays into the 
hidden sector while the curvaton decays into the standard model particles. Thus, the apparent 
fine-tuning of the initial condition of the NG curvaton based on the conventional naturalness 
argument may be due to our ignorance of the mechanism about how those parameters are 
determined. 

Another important constraint on the curvaton model is that it must account for the 
red-tilted power spectrum, ra^ — 1 ~ —10^^. As long as the curvaton potential is time- 
independent, this implies that the curvature of the curvaton potential is only about one order 
of magnitude smaller than the Hubble parameter during inflation. Then, unless the curvaton 
is initially sitting near the hilltop, it starts to oscillate about the potential minimum soon 
after inflation, and the curvaton mechanism does not work for most inflation models. Thus 
we are led to consider the hilltop curvaton models. However, the situation changes if we allow 
the curvaton potential to depend on time. In particular, the curvaton may have a negative 
Hubble mass during inflation, V D —nH^a^. Then the red-tilted spectrum simply implies 
K, ~ 10^^. After the infiaton decays, the negative Hubble mass term disappears, and the 
curvaton potential may become very flat, and the curvaton starts to oscillate long afterward. 
In this case we expect that the curvaton scenario works for many inflation models while 
accounting for n^ — 1 ~ — 10~^.^^ 

^^Wc note however that the curvaton dynamics would be comphcated with a negative Hubble mass term. 
For instance, the curvaton slow-rolls on the negative Hubble mass away from its origin during inflation, then 
long after reheating, the curvaton comes back and starts to oscillate about the origin. Alternatively, the red- 
tilted power spectrum can be easily realized in the curvaton model with a global U(l) symmetry; if the radial 
component of the complex scalar field moves from a large field value to a smaller one, the phase component 
(i.e. NG curvaton) would have a scale-dependent fluctuation, which results in the red-tilted power spectrum. 
MK and FT thank Shinta Kasuya for discussion on this issue. 
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Lastly we comment on a possible reason why the NG curvaton has a mass close to the 
Hubble parameter during inflation. If the curvaton potential is generated by non-perturbative 
dynamics, e.g. strong gauge interactions, the potential height is expected to be determined 
by the dynamical scale, and it is exponentially sensitive to the value of the gauge coupling at 
the UV scale. So, if there is a bias toward a larger value of the coupling constant at the UV 
scale, the curvaton mass is favored to be as large as possible. If the curvaton mass exceeds 
the Hubble parameter during inflation, it is stabilized at the potential minimum, and the 
curvaton mechanism does not work. Thus, the curvaton mass slightly smaller than ifinf may 
be selected in the curvaton landscape. Of course this is just a theoretical possibility, and it 
is hard to prove that such an argument is indeed the correct explanation. 

5 Conclusions 

In this paper, we have explored density perturbations from a curvaton with a generic po- 
tential, and shown that the curvaton's fleld fluctuations can lead to non-uniform onset of 
the curvaton oscillation 6 Hose- This sources additional contributions to the resulting density 
perturbations, in addition to that from the fleld fluctuations directly perturbing the curvaton 
energy density. We showed that the combination of the two effects gives rise to non-trivial 
behaviours for density perturbations from curvatons, such as the curvaton sourcing large /nl 
with either sign even if it dominates the universe before decay. The effect due to SHosc be- 
comes important when the curvaton potential deviates from the familiar quadratic type. Since 
a curvaton potential with negative curvature is required during inflation in order to explain 
the observationally suggested red-tilted density perturbation spectrum (without invoking in- 
flation with large \H/H'^\), it is of importance to evaluate effects from SHosc- Our work was 
also motivated by curvaton models based on microscopic physics such as supersymmetry or 
string theory, which provides explicit examples of curvaton potentials possessing non-trivial 
forms. 

As a simple case where the SHosc effects become dominant on the generated density per- 
turbations, we have studied curvatons that start their oscillations from a flat region along its 
potential. Especially for hilltop curvatons, we showed the strong enhancement of the linear 
order density perturbations Vq, accompanied by a mild increase of the non-Gaussianity /nl- 

Another non-trivial case, which we have not studied intensively in this paper, is when 
various contributions to the density perturbations cancel each other and thus suppress the 
perturbation amplitude. Conversely, then the non-Gaussianity is expected to become large. 
It would be interesting to investigate such cases, which may be realized by, e.g., curvaton 
potentials possessing plateau regions, or with superimposed repeated modulations.^^ We 
leave this for future work. 



^■^Superimposed repeated modulations to the inflaton potential can also have interesting consequences to the 
density perturbation spectrum generated by the inflaton itself, see e.g. [43, 44, 45, 46, 47, 48, 49, 50, 51, 52]. 
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We also applied our generic results to study a simple model where the curvaton is a 
pseudo-Nambu-Goldstone (NG) boson of a broken U(l) symmetry. Although NG curvatons 
away from the hilltop require high inflation/reheating scales in order to generate red-tilted 
density perturbation spectra with n^ — 1 ~ — 10~^, it was shown that NG curvatons at the 
hilltop can work with a wide range of inflation/reheating scales. We also showed that working 
NG curvatons in the hilltop predict the non-Gaussianity to lie in the range 10 < /nl ^ 30, 
which is accessible to upcoming CMB observations such as the PLANCK satellite. 

For the four-point correlation functions of the density perturbations from hilltop NG 
curvatons, one obtains relations that are similar to those for the familiar quadratic curvatons, 
i.e., tnl ~ 5'NL ~ /nl- However, different curvaton potentials may lead to different behaviours 
among the non-Gaussianities (see related discussions in [12]). It will be interesting to examine 
systematically the range of possibilities arising from curvatons with generic energy potentials. 

Both from observational and theoretical reasons, it is important to investigate in detail the 
conversion process of the curvaton field fluctuations into the resulting density perturbations. 
In this work, we have provided generic analytic expressions incorporating effects due to non- 
uniform onset of the curvaton oscillation, which open up new possibilities for the curvaton 
paradigm. 
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A Scalar Field Dynamics in an Expanding Universe 

In this appendix we discuss dynamics of a scalar field a rolling along a potential V{(t) in an 
expanding universe background. Here the potential V{a) is assumed to be a function only of 
0", and does not depend on say, time, explicitly. We especially show a condition under which 
the equation of motion of the field 



is approximated by 



where c is a constant. 



d + iH{t)d+?^^ = (A.l) 



cH(t)a . -^. (A.2) 
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The background expansion is considered to be independent of the cr-dynamics (as in the 
case of a curvaton before it starts to oscillate), and the universe to be dominated by some 
sort of energy density with equation of state p = wp (with constant w), so that the Hubble 
parameter satisfies 

^-^- (") 

By taking a time derivative of both sides of (A. 2) and comparing with (A.l), one can check 
that 

where a prime denotes a derivative with respect to a. Hence for 



V" 



cH^ 



< 1, (A.5) 



we obtain 

c^3+^ ,A.6) 

which is no less than 3 for w > —1. For example, c = 3 for a de Sitter, c = 9/2 for a matter 
dominated, and c = 5 for a radiation dominated universe. (We note that a non-zero but 
constant IVl/H"^ can also realize the approximation (A. 2) with different values of c, as in 
rapid-roll inflation [53, 54, 55, 56, 57, 58, 59].) 

One can further show that (A. 2) is an attractor. Introducing 



where we now define c by (A. 6), then one can check that 

c^ - c^^ (A.^ 



i y" .. ..2 



H cm 

Since c given by (A. 6) is positive, this equation shows that, at least when |^| < 1, |,^| damps 
as the universe expands until its amplitude becomes as small as C(|y'|/c^iJ^). Hence as long 
as the condition (A.5) holds, the approximation (A. 2) with (A. 6) is a stable attractor. 

B Density Perturbations from Curvatons with Non- 
Sinusoidal Oscillations 

In this appendix we generalize the expressions for density perturbations in the main body 
of the paper to incorporate a period of non-sinusoidal curvaton oscillation during which the 
energy density redshifts as pa oc a^" with arbitrary ra, followed by the usual sinusoidal os- 
cillation with n = 3. The results here can be applied to, e.g., a curvaton potential which is 
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quadratic around its (local) minimuni, but is dominated by a non-quadratic polynomial term 
at large a. 

Instead of the assumptions in Subsection 2.1, here we suppose that at some time tosc after 
inflation the curvaton starts an oscillation which redshifts away its energy density as p^ oc a~" 
with an arbitrary constant n, until its energy density becomes as small as pa = pasm, then the 
curvaton suddenly switches to a sinusoidal oscillation giving p„ oc a~^. After some period of 
the sinusoidal oscillation, the curvaton decays to radiation. Having in mind potentials that 
switch between quadratic and non-quadratic polynomials at a certain field value, we take the 
energy density pasm as a constant which is independent of a*. We also consider the curvaton 
energy density fraction to be negligibly small until treh or tosc, whichever is later. 

Then one can compute the SM' similarly as in Subsection 2.1, yielding 

dAf r (9 /3 \ 

Inp^osc-^lni/L, (B.l) 



da^ 4 + 3r da^ \ n 
d^Af 16(1 + r) fdAfV r d^ f 3 



+ -TTW- 1^ - In P<xosc -A\n Hl^ , (B.2) 



dal (4 + 3r)r \da-^ J 4 + 3r Scr^ \n 

where A = 3/4 for treh < ^osc, and A = 1 for tj-eh > tosc (in this latter case treh can be 
either before or after tgin, but of course before tdec)- We note that (2.20) holds in the present 
case as well. It is clear that the above expressions reproduce those in Subsection 2.1 for 
n = 3. Moreover, the discussions in Subsections 2.2 and 2.3 can be applied to the present 
case without any modification. 

example: Linear Curvatons 

As an example for applying the above formulae, let us take a look at a curvaton potential of 
the form 

2" 



n-) = AMf ' 



-1/2 

(B.3) 



with constant mass scales A and /. The potential is quadratic around its minimum a = 0, 
but approaches a linear potential at |o-| ^ /. 

Considering the curvaton to start its oscillation in the linear regime, i.e. cTosc ^ /, then 
by solving (2.31), one finds 

2c- 6 / AC A* \ ,„,, 

The energy density of the curvaton oscillating along the linear regime redshifts as p^ oc a~^, 
hence one can use the formulae (B.l) and (B.2) to obtain 



d^ (3 + 2A)r f K A^ 



da, 2(4 + 3r) V 3 fH,^, 
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-1 



^*-^TT^] ' (B-5) 



5 / 16(l + r) 2(4 + 3r) ) 

^^^"e^^l^Ts;^ 3T2:4"/' ^^-^^ 

which exhibit behaviours that are rather close to those of curvatons with quadratic potentials. 
However, here the spectral index (2.45) vanishes at leading order, when H^ = 0. 
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